Dark stationary matter-waves via parity-selective evaporation 
in a Tonks-Girardeau gas 
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We propose a scheme for observing dark stationary waves in a Tonks-Girardeau (TG) gas. The 
scheme is based on parity-selective dynamical "evaporation" of the gas via a time-dependent po- 
tential, which excites the gas from its ground state towards a desired specially-tailored many-body 
state. These excitations of the TG gas are analogous to linear partially coherent nondiffracting 
beams in optics, as evident from the mapping between the quantum dynamics of the TG gas and 
the propagation of incoherent light in one-dimensional linear photonic structures. 
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Trapped bose gases confined to one-dimensional (ID) 
geometry are highly attractive for studying quantum 
many-body dynamics because pertinent models yield ex- 
act solutions 0000; these regimes are experimentally 
accessible 0,0013; and quantum effects are enhanced 
Ell EH]- Tonks-Girardeau (TG) gas is a system of 
ID bosons with "impenetrable core" repulsive interac- 
tions lj. In the fashion of the Pauli exclusion princi- 
ple, "impenetrable cores" prevent bosons to occupy the 
same position in space, which causes TG gas to exhibit 
fcrmionic properties. This similarity is manifested in the 
mapping between ID noninteracting fermions and TG 
bosons 0, 0. From the properties of atomic interac- 
tions in tight atomic waveguides ;9] it follows that the 
TG regime can be reached at low temperatures, low lin- 
ear densities or stronger effective interactions |9j, 1 1 0L l l l| . 
The experimental realizations of the TG gas in 2004 [aH 
boosted the physical relevance of the model. A recent ex- 
periment has demonstrated that such gas does not relax 
to the thermodynamic equilibrium even after numerous 
collisions |8(, due to the integrability of the underlying 
model. This work inspired a theoretical study of ID 
" impenetrable core" bosons on a lattice, which suggested 
that the system can undergo irreversible relaxation to a 
steady state carrying more memory of the initial con- 
ditions than the usual thermodynamic equilibrium |l2| . 
Some aspects of the TG quantum dynamics have been 
studied theoretically within the context of so-called " dark 
solitons" 

E3 

, matter- wave interference 

Q; 

ID expan- 
sion [ThL llfij l , irregular dynamics , and coherent states 
|l8| . In this work, we study the dynamical tailoring of 
the TG gas via a time-dependent potential to produce 
dark stationary states, as well as point out the relation 
between TG dynamics and the propagation of incoherent 
light in linear photonic structures. 

The Fermi-Bose mapping 0, , applicable both in the 
static and time-dependent case , prescribes the con- 
struction of the exact many-body wavefunction of the 
TG gas from single-particle (SP) wavefunctions, which 



obey a set of uncoupled linear SP Schrodinger equations. 
In Ref. 0] Girardeau and Wright discuss the dynam- 
ics of the TG gas within the context of dark solitons 
Dark solitons are fundamental nonlinear 
excitations, which have been studied theoretically within 
the nonlinear mean-field theories applicable for weakly 
interacting gases |l9l I21I l22l , and observed experimen- 
tally in these regimes 20]. For a strongly interacting 
TG gas on a ring, Girardeau and Wright noticed 
that if the many-body wavefunction is constructed solely 
from the odd-parity SP eigenstates of the system, the 
SP density will have a dip at zero, similar in structure 
to dark-solitons. However, such a specially structured 
many-body state is unlikely to occur without deliberate 
preparation, since even and odd parity SP states of that 
system are intermingled when ordered with respect to 
energy (see the discussion in Ref. 0). In the study of 
Busch and Huyet 0] , the collapses and reappearances of 
TG dark solitonlike structures in an harmonic trap are 
attributed to the mixture of the odd and the even compo- 
nents in the excitation. Generally, the SP eigenstates in 
parity-invariant ID potentials can be chosen to be either 
even or odd, which makes them candidates for observing 
dark stationary structures in the TG gas. However, for 
their experimental realization under such confinement, it 
is essential to separate components of different parity. 

Here we propose a scheme for observing dark station- 
ary waves in a TG gas. A time-dependent potential is 
used to selectively "evaporate" the even component of 
the many-body wavefunction, thereby creating a dark 
stationary wave. Such excitation of the strongly inter- 
acting TG gas is in fact an excited many-body eigenstate 
of the system, which distinguishes it from dark-solitons 
of the nonlinear mean-field equations applicable for weak 
interactions 0, [2£|, 0, • We point out that such ex- 
cited eigenstates of the TG gas are analogous to linear 
partially coherent nondiffracting beams in optics [^E^ . 
as evident from the mapping between the quantum dy- 
namics of the TG gas and the propagation of incoherent 
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light in one-dimensional linear photonic structures, pre- 
sented in this Letter. 

We consider N impenetrable bosons, confined within 
a ID external potential V ext (x,t). The fully sym- 
metrized many-body wavefunction describing the system, 
iPb(x\, . . . , xnJ^jIs constructed according to the Fermi- 
Bose mapping UQ. Let ?/>(£, r) denote a set of orthonor- 
mal SP wavefunctions obeying the set of uncoupled linear 
Schrodinger equations, 



.dipm 



ip m (€,r), m = 1, 



,N. 



(1) 

In order to unify notation and discuss the equivalence 
with optics, we find it convenient to use dimensionless 
units. The boson mass m and the (arbitrary) choice of 
spatial lengthscale xq (£ = x/xq) determine the units 
of time to = 2mx 2 ,/?i (r = t/to) and energy eo = 
h 2 /(2mxl) [V(£,t) = V e xt{x,t)/e }. From the SP wave- 
functions ip m one first constructs a fully antisymmetric 
(fermionic) wavefunction in the form of the Slater deter- 
minant, ip F (xx,...,XN,t) = \Jx$ /N\ det [tp m (£j , r)] ; ip F 
describes a system of spinless noninteracting fermions in 
the ID potential V ext (x, t) 0,0. The bosonic many-body 
solution ipB {%it ■ ■ ■■> xn, t) is obtained after symmetriza- 
tion of ipp: 



N N 



iJjb = A(xi,...,x N )\ -g- det M m (£j,T)}, (2) 

V iv m,j — 1 

where A = ni<i<-j<7vsgn(xi — x.j) is a "unit antisym- 
metric function" [lj. Thus, the quantum dynamics of 
the TG gas is obtained from Eq. J2Jl after solving Eq. 
Q . For example, the evolution of the single-particle den- 
sity p S p(x,t) = J dx 2 . ■ .dx N \ip B (x,X2, ■ ■ .,x N ,t)\ 2 cor- 
responds to the evolution of t) = ^2 m \ip m (i, t)| 2 
It should be noted that while some quantities of the cor- 
responding fermionic system are identical (e.g., the SP 
density), some significantly differ (e.g., the momentum 
distribution 0,0)- 

Excited many-body eigenstates with soliton-like SP 
density are found in real, time-independent, and par- 
ity invariant potentials, V(£) = V(— £). Let </> e , 7 (£) de- 
note the eigenstates, and let e denote eigenvalues (en- 
ergies) of the SP Hamiltonian H = -d 2 /d£ 2 + V(£) 
with appropriate boundary conditions. The extra in- 
dex 7 is used in case there are degenerate SP eigen- 
states. Since the Hamiltonian commutes with the par- 
ity operator, they can have a common complete set of 
eigenstates, in which case the eigenstates </> e , 7 (£) are ei- 
ther symmetric </>+ 7 (£) = <j>t,y(~Q (even parity) or an- 
tisymmetric <ftZ-y(£,) — ~^<T 7 (~£) (odd parity). Con- 
sider a bosonic many-body wavefunction, constructed ac- 
cording to Eq. ©, from odd-parity eigenmodes only, 
i.e., every ip m equals to one of the eigenstates <f>J^{£.) 



(V-Vn 7^ ipn f° r m =/= n). Such a wavefunction is an ex- 
cited many-body eigenstate of the system. Its SP density 
p~(£) = YH^7,-y(0\ 2 1S stationary, and p~(0) = be- 
cause 07 7 (O) = 0. The SP density of this excited many- 
body eigenstate thus has a dip at £ = 0, which resembles 
the structure of nonlinear dark-solitons i3, 110. 
Therefore we will refer to these states as dark many-body 
eigenstates. If we construct waves ?/> TO from symmetric 
eigenmodes </>+ 7 (£) only, the density p+(£) = J] |^ 7 (£)| 2 
is likely to have a pronounced peak at £ = 0, reminiscent 
of anti-dark solitons. Hence, they will be referred to as 
anti-dark many-body eigenstates. 

Dark and anti-dark many-body eigenstates are spe- 
cially tailored excitations of the TG gas, constructed 
solely from the odd or even SP eigenstates of the par- 
ity invariant potential V(t;). Assuming the absence of 
degeneracy, the even and odd parity SP eigenstates al- 
ternate when ordered with respect to energy, meaning 
that such specially tailored states are unlikely to natu- 
rally occur. We illustrate this by studying N = 20 TG 
bosons in the external container-like potential V c (£) = 
V r c {2 + E i = 1 ,2(-) i+1 tanh^(H(-)^)} (V c ° = 15, 
x w = 4, and x c = 7), shown in Fig. da). The SP 
density p of the ground state is plotted as a solid black 
line, whereas the dotted red (dashed blue) line depicts 
the odd (even, respectively) component of the SP density: 
p{Q = P + (£) + P~ (0 ■ The energies of the even (odd) par- 
ity SP eigenmodes are shown as blue squares (red circles, 
respectively). As expected odd and even parity eigen- 
states alternate with increasing energy. Thus, half of the 
SP eigenstates comprising the ground state are even and 
half are odd. The even p + (odd p~) component of the 
SP density has the structure of the anti-dark (dark, re- 
spectively) many-body eigenstates. In what follows we 
propose a method for the dynamical excitation of dark 
stationary waves. 

Our scheme utilizes a time-dependent potential which 
tailors the many-body wavefunction in a specific desired 
fashion, and separates the odd from the even compo- 
nent. The system of N = 20 TG bosons is initially in 
the ground state of the confining potential, see Fig^a). 
In the spirit of Refs. 0, , we perturb this system at 
its symmetry point £ = 0, with a spatially-narrow re- 
pulsive potential, which may be obtained with a laser 
0,0. However, in contrast to 0,0, here we period- 
ically switch this potential on and off in time. Such a 
time-dependent potential can be modeled as V'(£, r) = 
V r '[sgn(sin(27ri/r p )) + 1)] exp[-(£/er) 2 ], where r p = 0.1 
is the periodicity of the laser signal, Vq = 100 corre- 
sponds to the peak intensity of the laser, and a = 0.06 to 
its spatial focusing width. As bosons are kicked by the 
time-dependent potential, they acquire energy which has 
a certain probability of being higher than the lip of the 
trap and are ejected away from it. However, because the 
laser is focused close to £ = 0, it strongly affects only the 
even-parity SP eigenstates, whereas the odd-parity eigen- 
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FIG. 1: (color online) Dark and anti-dark many-body eigen- 
states: structure and excitation, (a) The container potential 
V c (£) (black dot-dashed line), the even p + (blue dashed line), 
odd p~ (red dotted line), and total SP density p = p + + p~ 
(black solid line) of the ground state with 20 bosons; the first 
20 eigenvalues are shown as blue squares (even SP eigenstates) 
and red circles (odd SP eigenstates). (b) The spectrum Pi of 
the excitation at r = and r = 19.8, and (c) the dynamics 
of the SP density following the time-dependent perturbation, 
(d) The SP density p(£) of an anti-dark (black dotted line) 
and dark (red solid line) many-body eigenstate in the periodic 
potential (see text for details). 



states are left nearly unperturbed. Consequently, the 
many-body wavefunction within the container < x c ) 
takes on a specific structure: it is constructed via Eq. 
mainly from the odd-parity SP eigenstates. This filtering 
process is depicted in Fig. ^b), showing the spectrum 
of the SP wavefunctions ip m (£,, t) calculated according to 

= E m I JX d$i> m (Z,T)<t>t(0\ 2 , where is the 
ith eigenstate of the SP Hamiltonian. The spectrum at 
t = is flat (black squares) because the odd and the even 
eigenstates are equally present. However, the spectrum 
after r = 19.8 (red circles) is mainly comprised from 
odd-parity eigenstates. Fig. [He) shows the evolution of 
the total SP density. The time-dependent potential acts 
within the interval r = [0, 19.8]. After r = 19.8 (marked 
by a horizontal line) it is turned off. The SP density 
nevertheless retains a dark notch at £ = even after the 
time-dependent potential is turned off, clearly displaying 
dark stationary wave evolution. The numerical evolu- 
tion of Eq. (JTJ is performed with the split-step Fourier 
method. 

Although different in nature, the scheme proposed here 
is akin to evaporative cooling. The concept should work 
for various types of container-like potentials. The scheme 
is also fairly robust. For lasers with larger intensity oc Vq, 
the filtering occurs on faster time scales (i.e., a smaller 
number of on-off switches are sufficient) and is more effi- 
cient. The focusing width of the laser a limits the number 



of particles which can be efficiently filtered. This width 
should be sufficiently smaller than the period of the spa- 
tial oscillation (close to £ = 0) of the Nth SP eigenstate. 
It should be emphasized that even though we follow the 
spirit of Refs. [j,[l3|j the proposed scheme separates the 
odd and the even component in space, while the many- 
body wavefunction within the container assumes the par- 
ticular structure of a dark stationary wave. 

While the proposed method for exciting dark station- 
ary states of the TG gas employs a container- like po- 
tential, it should be emphasized that the notion of dark 
and anti-dark many-body eigenstates pertains to various 
parity invariant potentials. We illustrate this fact in a pe- 
riodic potential V(£) = V(£ + D) (e.g., optically induced 
lattices). The SP eigenstates of this system are Bloch 
waves [25j of the form 0a i71 (£,t) = Ufc,„(£)e lfc? e" I£fc '" T , 
where n denotes the band number, k is the Bloch wave- 
vector, and Uk, n (€) = u k,n(£, + D) describes the periodic 
spatial profile of the Bloch wave. Since V(£) = V(— £), 
the Bloch waves <pk,n and 4>-k, n are degenerate. By prop- 
erly choosing the coefficients a km and a_fc m within the su- 
perposition ip m (£,T) = a km 4) kmtn + a-k m (f>-k m ,m degen- 
erate eigenstates <^±fc m ,„ can be superimposed to obtain 
even [4>+(£,t) = ^+(-£,r)] and odd-parity [ip~(£,T) = 
— ip m (— £> T )] eigenstates. The many-body wavefunction 
comprised solely from r) [ip m (^, r)] via Eq. J2J is a 

dark (anti-dark) excited many-body eigenstate of the TG 
gas in the lattice. Figure ^Jd) shows the (stationary) SP 
density of the dark and anti-dark many-body eigenstate 
in a periodic potential V(£) = 10cos 2 (7rx), constructed 
by symmetrizing the lowest N = 21 SP eigenmodes. The 
calculation is performed on the ring (periodic boundary 
conditions) of length L = 71. 

Relation to linear incoherent light. - Dark and anti-dark 
excited many-body eigenstates of the TG gas are analo- 
gous to partially coherent nondiffracting beams that were 
studied in the context of classical optics |23L 0|. In 
order clarify this point, we first demonstrate the map- 
ping between the propagation of incoherent light in lin- 
ear ID photonic structures and the TG gas dy- 
namics. Consider a quasimonochromatic, linearly po- 
larized, partially-spatially incoherent light beam which 
propagates paraxially in the ID photonic structure de- 
scribed by the spatially dependent index of refraction 
n tot — n,Q+2non(x, z). The classical electromagnetic field 
E(x, z, t) of the beam randomly fluctuates; z [x] denotes 
the propagation axis (spatial, respectively) coordinate. 
The state of the system is described by the mutual co- 
herence function B = (E*(x 2l z, t)E{x\, z, t)) H3, where 
brackets denote the time-average, which equals ensem- 
ble average assuming the light source is stationary and 
ergodic |27| . The mutual coherence function B can be de- 
composed through an orthonormal set of coherent modes 
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ip m and their modal weights A m |27 



B(xi,x 2 ,t) A m V40r 2 , z)ip m (xi, z). 



(3) 



in optics. This analogy is a consequence of the mapping 
between incoherent light in linear ID photonic structures 
and the TG gas. 



In order to connect B to Eq. we switch to dimen- 
sionless units: £ = x/x , r = z/(2kXo) ("time" is here 
the propagation length) where k = tiqui/c, and u> is the 
temporal frequency of the beam. The potential V arises 
from the refractive index V(£, r) = — 2{kxo) 2 n{x, z)/tiq. 
Waves_0 m (£,r) = y/x^i) m {x,z) obey Eq. JJ) (e.g., see 
Ref. |2q). affirming the mapping between the two sys- 
tems. Note that each solution of Eq. 0) generates one 
bosonic many-body wavefunction via Eq. J2J), and, due 
to the arbitrary choice of modal weights A m , many corre- 
lation functions (J3J corresponding to incoherent optical 
fields propagating in linear ID photonic structures. The 
density p = J2 m IV'ml 2 corresponds to time-averaged in- 
tensity / = J2 m AmlVU 2 E3 ■ 

One particular example of a partially coherent non- 
diffracting optical beam propagating in vacuum, cor- 
responds to the dark stationary TG wave on a ring 
studied by Girardeau and Wright 0. An incoher- 
ent optical beam with the mutual coherence function 
B(xi,X2) = J dk x G(k x ) sm(k x x), or equivalently with 
the modal structure ij)k x {%) — y G(k x ) sm(k x x), is propa- 
gation invariant. If its power spectrum G(k x ) is rectangu- 
lar, G(k x ) = Iq/K for \k x \ < K, and zero otherwise, the 
intensity structure has the form Iq/2[1 — jo(Kx)], which 
is exactly the form of the odd SP density-component of 
the dark stationary wave on a ring in the thermodynamic 
limit Q. If B(x\, X2) = J dk x G{k x ) cos(k x x), one obtains 
anti-dark optical propagation-invariant waves. 

We have thus established a link between the dynamics 
of incoherent light in linear photonic media and TG gas 
via Eq. JIJ. It should be kept in mind the the former 
system is classical, while the latter is quantum, and the 
evolution of the quantities derived from the set of waves 
ip m (e.g., the density p = J2 m IV'ml 2 ) should be properly 
interpreted. This mapping adds to the analogies between 
optical and matter waves |29j |. and in particular to the 
analogy between nonlinear partially coherent optical- and 
matter- waves |.3C)j. 

Before closing, it should also be noted that the 
evolution equation for the mutual coherence function 
B(x±, X2,z), in the paraxial approximation, in linear pho- 
tonic structures (e.g., see |28|), is identical to the evolu- 
tion of the reduced single-particle density matrix of non- 
interacting spinless fcrmions (in ID, and 2D as well). 

In conclusion, we have proposed a scheme for excit- 
ing dark stationary waves of the TG gas. Within our 
scheme, a time-dependent potential focused to the center 
of the trap, selectively "evaporates" a non-desirable part 
of the many-body wavefunction, thereby creating a dark 
stationary wave. The stationary waves of the TG gas 
are analogous to partially-coherent nondiffracting beams 
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